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SOLUTION OF PROBLEM 155. 



BY CHAS. H. KUMMELI., V. S. LAKE SUEVEY, DETROIT, MICH. 

"To find the least distance between two places given by latitude and lon- 
gitude, taking into account the polar compression." 

This question admits two interpretations, viz.: — 1, To find the rectilinear 
distance, that is the chord joining the two points; and 2, To find the short- 
est distance between the two points on the spheroid or that along the geo- 
desic line. 

First 8xi/pposition. — Let cp^, tp^ be the latitudes and X-^, l^ ^^ longitudes 
of the two given points respectively then we have their orthogonal coordi- 
nates [x-^, 2/i, z^ and {x^, y^, z,^) as follows: 

Xy^ = iV^jCos f jCosAj, Xc^ = iVjCos ^2^03^2/ 

2/1 = iV^iCosf icos/ii, 2/2 = N^COBf^COSi.^, ^ (1) 

«! = iVi(l-e^)sin9?5L, «2 = -^2(l~^^)sinf2 






where N = — 7- tt-^^-—^ = normal ending at polar axis, a = equatorial 

1/(1 — e^sm^) ox- ) 1 

semiaxis and e = eccentricity of meridian. We have then, denoting the 
straight line connecting the two points by s, 

s = l/[(a;2— «'i)'+(2/2— 2/i)'+(%— «l)'] 
= \/\_N\ cos^f 1 — 2N.^N^cos f ^ cos ^j cos (^2 — •^1) + -^1 cos'f 2 

■^{l—e')\N^sm <p^ -N^sin <p^f] 
= ]/-( iVf — 2NiN2[cos<piCos<p2Cos{?.2 — ^Ai)-|-sin^isinyi2] 

+Nl-^{2-e'){N^sm ,p^—N^sm (pj )■ (2) 
If we draw through the center of the spheroid two lines parallel to the 
normals i\^i, N^ we have, denoting their included angle by a, 

cos a = sin ^^^sin ^2 +cos ^^ cos f^ cos (^2 — ^i)- (3) 

We have then 
s = ylNl—INT^N^cfB ff+Nl—e^i^—^^N^sm f^—N^sin <p{f]. (4) 
Join the ends of the normals by a straight line which may be denoted by s', 
then we have from the rectilinear triangle thus formed 

s'2 = Nl—2NiN^cos a+Nl (5) 

hence 

s = y'ls'^-e\2—e'){N^sin f^—N^sin <p^)'^. (6) 

The geometrical signification of the second term under the radical, which is 
always a small quantity and may be neglected for two places of nearly the 
same latitude, is the following : 



—US- 
Let Az' =■ N^ws. <p^ — iVisin ^1 = projection of s' on polar axis, 

Az = ^2 — ^\ = (1 — fi^)(-^2si^f 2 — -A^isin^j) = projection of s on 
polar axis, then Az'^ — Az^ = e\2 — e^)(iV2sin ^^ — -^isin (p-^J and 

s = i/ls'^—Az'2+Az'l (7) 

hence s'^—A-/ = s'^—Az'"^, (8) 

that is, the projections of s and s' on the equator are equal. 

Second Supposition— Any surface of revolution may be represented by the 
equation z =f{p), (1) 

where z = distance of any point from some plane of reference at right an- 
gles to the axis of rotation, which may be called the equator, and p = radi- 
us of parallel. In this form the equation of the surface is independent of 
the third coordinate, the longitude ^ = inclination of any meridian to the 
first meridian. In this system of coordinates we have for the element of arc 
of any curve in space 

ds = ^ldp'+pm^+dz% (2) 

and the arc of any curve between two fixed points {p^, X^, Zj) {p^, X^, z^) 

s = f^-i/ldp'+pm^-\-d^l (3) 

If this arc is to be on the surface z = f{p) there will be only two inde- 
pendent variables, viz., p and A in (3) ; and if we introduce the condition 
that s shall be the shortest arc on the surface between the given points we 

must have -5^ = 0, or else -5- = 0. 
0/ op 

The first term vanishes at the limits since the curve is to pass through fixed 
points, consequently the other term must be also, and therefore 

d J = 0, or, integrating p'dk = eds, (4) 

and integrating again 

PpHX = cs. (6) 

This equation shows that the arc of the curve is proportional to the sector 
of its horizontal projection. 

If A is the angle which the geodesic line makes with the meridian or its 
azimuth, we have evidently 

sin^=/,| (6) 

Combining this equation with (4) we have 

^sinJ. = c. (7) 
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The constant c is therefore the radius of parallel where the geodesic line 
meets a meridian at right angles. Combining (4) with (2) we obtain 

which is the polar differential equation of the horizontal projection of the 
geodesic line on any surface of revolution ; for it is an equation between the 
longitude )., the radius of parallel p and an undetermined constant c. To 
determine this constant we have the condition 

The length of the arc of the geodesic is then 

For the spheroid we have 

z = vK^-e^{o?-p')'\; (11) 

The latitudes and longitudes of the two points being the data, we have for 
the limits of these integrals 

jOj = a(cosf i)-^if 1, |02 = a (cosp2)-^'^^2) (12) 

if A(p =z -[/(l — e^sinY). 

Placing p = |/[c^+(a'' — c^)sinY] we have 

^ ^ -^ 1 o^+{d' — c^)smY ^/(a^— cVj ] L 

s = /V[«='-cV-(a^-c^)e%inV]=l/(a»-cV)^' [^, e^(^£^^)] 

' . • • • (^^") 

These elliptics give the complete solution of the problem in theory. The 

great difficulty consists in the determination of the constant c. If the azi- 
muth of the geodesic were identical with the azimuth as measured by a 
theodolite and which is the inclination of the vertical circle, passing through 
the other point, to the meridian;* then equation (7) would give the constant c 

*This field azimuth is given in terms of the latitudes and longitudes of the two points by 

the formula: cot^i =cot a, - ff^^'t i.Mt. \^J^-^^] , where cot a^ = [tan f, 

sin(V^i)cos(i2 I A<»2 Ai-i J 
Xcos ^j-sin ^jcos (/l2-/li)]-s-sin(/l2-/'.j). It may differ from the azimuth of the geodesic line 
by any amount between 0° and 90°. For instance, if both points are on the equator and 
nearly 180° apart the geodesic line passes near the pole, its azimuth is therefore nearly 180° 
while the field azimuth (of course impr'le) would follow the plane of the equator and=270°. 
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without any difficulty and equation (9) would be superfluous. If the two 
points are not too distant a value for o from (7) may be used as a first ap- 
proximation and its correction determined by (9). The constant o being 
thus found, the length of the arc is found by the evaluation of an elliptic 
of the second species. 

Note. — As it is contended that in the published answer to Prof. Hall's 
Query (see p. 94) the series which represents the value of u converges so 
slowly that the method is inconvenient, another answer is here submitted 
as given by Chas. H. Kummell. 

To find the most convenient way of computing the numerical value of the 
definite integral 

1= f^dw-i/ismip). (1) 

We have 

/fd^(sin^)2"-i(cos^)2''-i = i5(m,n) = ff^^y (2) 

Placing m =: f and w = f we have 

1= pd<py^{sm<p) = ^fjf^^f . (3) 

But r(j) = i/tt 

r(|) =-^ by the theorem: r{n)r{l—n) = 



therefore 



r{l) • \ J \ sin n;r 

r(|)=ir(i)"" " : r{n+l) =nr{n); 

T- (M 



The circumference of the lemniscate of Bernoulli may be expressed in 
terms of r{l) as follows: 

The polar equation of the lemniscate is r^ = a' cos 2^. Denoting the 
circumference by P we have 

Placing 6 = Itt — |^ we have 



2^„T/(sin^) * '"'' 4- rH) 
P = 
Combining this with (4) we obtain 



^=«7W' '") 



J= ?H. (7) 



